Method for direct observation of coherent quantum oscillations in a superconducting 

phase qubit. Computer simulations. 
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Time-domain observations of coherent oscillations between quantum states in mesoscopic super- 
conducting systems have so far been restricted to restoring the time-dependent probability distri- 
bution from the readout statistics. We propose a method for direct observation of Rabi oscillations 
in a phase qubit. The external source, typically in GHz range, induces transitions between the 
qubit levels. The resulting Rabi oscillations of supercurrent in the qubit loop are detected by a high 
quality resonant tank circuit, inductively coupled to the phase qubit. Here we present the results of 
detailed computer simulations of the interaction of a classical object (resonant tank circuit) with a 
quantum object (phase qubit). We explicitly account for the back action of a tank circuit and for 
the unpredictable nature of outcome of a single measurement. According to the results of our sim- 
ulations the Rabi oscillations in MHz range can be detected using conventional NMR pulse Fourier 
technique. 

PACS numbers: 03.65.Ta, 73.23.Ra 



I. INTRODUCTION 



As is known the persistent current qubit (phase qubit) 
is one of the candidates as a key element of a scalable solid 
state quantum processor The basic dynamic manifes- 
tations of a quantum nature of the qubit are macroscopic 
quantum coherent (MQC) oscillations (Rabi oscillations) 
between its two basis states, which are differed by the 
direction of macroscopic current in the qubit loop. 

Up to now Rabi oscillations in the time domain2i4 or as 
a function of the perturbation power^ have been detected 
indirectly through the statistics of switching events (e.g., 
escapes into continuum). In either case the probability 
P{t) or P{E) was obtained and analyzed to detect the 
oscillations. 

More attractive is a direct detection of MQC oscil- 
lations through a weak continuous measurement of a 
classical variable, which would implicitly incorporate the 
statistics of quantum switching events, not destroying the 
quantum coherence of the qubit at the same time&^ 

In this paper we describe an approach which allows a 
direct detection of MQC oscillations of macroscopic cur- 
rent flowing in a loop of a phase qubit. To be specific, we 
will use the example of three- Josephson-junction (3JJ) 
small-inductance phase qubit (persistent current qubit^) 
where level anticrossing was already observed»^ 

In our method a resonant tank circuit with known in- 
ductance Lt, capacitance Ct, and quality Qt is coupled 
with a target Josephson circuit through the mutual in- 
ductance M (Fig. III. 

The method was successfully applied to a 3JJ qubit 
in classical regime^" when the hysteretic dependence of 
qubit energy on the external magnetic flux was restored 
in accordance to the predictions of Ref. 0- 

Here we extend the approach which is described in our 
earlier paper. We explicitly account for the back ac- 
tion of a tank circuit and for the unpredictable quantum 
mechanical nature of outcome of a single measurement. 




FIG. 1: Phase qubit coupled to a tank circuit. 



According to the results of our simulations the Rabi oscil- 
lations in MHz range can be detected using conventional 
NMR pulse Fourier technique. 



II. QUANTUM DYNAMICS OF 3JJ FLUX 
QUBIT 

Quantum dynamics of this qubit has been studied in 
detail in Ref. ^ The qubit consists of a loop with three 
Josephson junctions. The loop has very small induc- 
tance, typically in the pH range. It insures effective 
decoupling of qubit from external environment. Two 
Josephson junctions have equal critical current /c and 
capacitance C, while the critical current and capaci- 
tance of a third junction is a little bit smaller, ale, aC, 
where 0.5 < a < 1. If the Josephson coupling energy 
E} = /c$o/27r, where $o = h/2e is a flux quantum, is 
much more than the Coulomb energy Ec = e'^/2C, then 
the phase of a Cooper pair wave function is well defined. 
As was shown in Refs. in the vicinity of $ = $o/2 
this system has two quantum stable states with persis- 
tent circulating current of opposite sign. Therefore, the 
persistent current qubit can be described by following 
two-level Hamiltonian: 



where tr^ and ax are Pauli spin operators. 



(1) 
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Hamiltonian Hq is written in the flux basis: the basis 
of localized (left, right) states, so that two eigenvectors 
of (Jz correspond to the two classical states with a left 
or a right circulating current. A is the tunnel frequency 
and £{t) is in general a time dependent bias which is 
controlled by externally applied flux 



$(t) =$, + $ac(t), 



(2) 



where $a; is a time independent external flux, '^ac{t) is a 
monochromatic high frequency signal from the external 
source. According to |2Jl we write e{t) in the form: 

e{t)=£a+e{t). (3) 

In the absence of time dependent flux {£{t) — 0) the 
eigenstates of Hamiltonian Hq are 



Notice that at the degeneracy point (/ = 0) the low 
frequency part of the average current, which is given by 
first term in vanishes. 

Below we define the density matrix elements poo = 
|C_(i)P, Pu = \C+{t)\\ pio = C+{t)C*_{t), poi = Pto> 
Poo + Pii = 1 which are related to spin operators: 
{■^{t)\az\^{t)) = 2poo - 1, {■^{t)\ax\'^{t)) = 2Repio, 
{^{t)\aY\^{t)) = 21m/9io. We take the high frequency 
excitation in the form e{t) = ei cos(wt) and rewrite 
Hamiltonian Q in the eigenstate basis '^±: 



Hq = -h in 



cos{Lut)j az+h^-^ cos{Lut)ax- (8) 



From density matrix equation ihp(t) = [Hq,p{t)] we 
get the following set of equations for the matrix elements 
of p: 



(4) 



dA Aei , 

— = 2—^ cos(ujt)B. 



where eo = EjfX{a,g)/n, f = $,/$o - 5 « 1- The 
explicit dependence of A(a, g) on qubit parameters a and 
g = Ej / Ec has been found in Ref.[l2- The fiux states are 

the eigenstates of a macroscopic current which circulates = — cos{ujt){2A — 1) 

in a qubit loop. Within a harmonic approximationi^ it 
is not difficult to find current operator in the fiux basis: 
Iq = Ic{\{a,g)/27r)az. 
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Hamiltonian Hq can be written in eigenstate basis as 
Hq = —Tiilaz with eigenfunctions 5'±: Hq^± = E±'^±. 
The stationary state wave functions ^± can be written 
as the superpositions of the wave functions in the fiux 
basis '^L,'^R, where L, R stand for the left and right 
well respectively: ± = a±'^ l + b±^ r, 



A 



:,6± = 



(5) 



The transformation from fiux basis to eigenstate basis 
is performed with the aid of rotation matrix R{C) — 
exp(icry^/2), where cos^ = £o/i^, sin^ = A/il. There- 
fore, the current operator in eigenstate basis is 



In 



\{a,g) 



(6) 



The high frequency excitation applied to the qubit in- 
duces the transitions between two levels which result in 
a superposition state for the wave function of the sys- 
tem: *(i) = C+(t)*+-|-C_(0*-. The coefficients C±(t) 
can be obtained from the solution of time dependent 
Schrodinger equation with proper initial conditions for 
C±{t). From © we obtain the average current in the 
superposition state ^'(i): 



dC 
H 



2 n 



cos(ti>t) ) B 



(9) 



cos(a't)j C, 
(10) 

(11) 



where A ~ poo, R — Impio, C ~ Rcpip. 

These equations cannot be solved analytically, how- 
ever, we can try to predict the evolution of the quantities 
A, B, C since the Eqs. (I9I1UI11I) are similar to those of 
a free spin in external magnetic field in NMR, where A 
is analogue of longitudinal magnetization Mz, and B, C 
are analogues of transverse magnetizations My and Mx , 
respectively. The only difference is that here the only de- 
tectable quantity is az component which is proportional 
to the circulating current. Therefore, we could expect 
the evolution of quantities A, B, C under the influence 
of external high frequency excitation is similar to that in 
NMR. If the frequency of external excitation u! is close 
to the gap frequency 2Q, the main harmonic of A will be 
the Rabi frequency 



LOR ^ ^{uj-2n)^ + n% 



(12) 



where flR depends on qubit parameters and on the am- 
plitude of excitation £1, while the quantities B, C will 
oscillate with the gap frequency 257, modulated by the 
Rabi frequency: B,C ^ sin ri/ji sin 2rit. If £1 <C 51 we 
can estimate ftji from Eqs. H9I10I11() in rotating wave 
approximation. We obtain 



/c^^(|C_(t)| 



-Ic 



2tt n 
A(a,g) A 
TT n 



|C+(t)P) 
Re{C+{t)C*_{t)). (7) 



A£i/57. 



(13) 



Since we want here to treat the problem exactly we solved 
Eqs. (|9I1UI11I) numerically. 
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FIG. 2: Time evolution of A and B for qubit without dissi- 
pation. 




t(^S) 



FIG. 3: Time evolution of A and B computed from Eqs. (14)- 
(16). Tr = 50 ms, Td = 200 ns. 



For all computer simulations we take the following pa- 
rameters of qubit: Ic = 400 nA, L — 15 pH, a — 0.8, 
g = 100, n/2Tr = 250 MHz, Eq = A. The external high 
frequency is equal to the gap: w = 2fl. The excitation 
amplitude ei = 0.28V. which ensures the Rabi frequency 
Q.ii/2ti = 50 MHz. In addition, we assume the qubit is 
in its ground state in initial time ti = 0, so that j4(0) = 1, 
5(0)=C(0)=0. In all simulations the time span is 10 /iS 
from which the particular time windows had been choose 
for the figures. 

The time evolution of A{t) and B{i) is shown on Fig.|21 

It is clearly seen that B oscillates with gap frequency, 
while the frequency of A is almost ten times smaller: 
(oscillation period of B: Tb ~ 2 x 10"^ s, while the same 
quantity for ^4 is Ta « 2 x 10~^ s. The small distortions 
on A curve are due to a strong deviation of excitation 
signal from transverse rotating wave form, while B curve 
is clearly modulated with Rabi frequency Q.^ (Fig.|21). 

The crucial requirement for the proper operation of a 
qubit is the preservation of phase coherence under influ- 
ence of dissipative environment. We include the environ- 
ment effects phenomenologically in Eqs. H9I10I11|I : 



^ = 2^ cos{ojt)B -Y^{A- Ao), (14) 



— = -^cos{ut)(2A-l)-2 



(15) 



cos{Lot) 



(16) 



where is the equilibrium value of density matrix 



^0 = Poo = 



1 + tanh 



no. 



Tr and Ta are relaxation 



and dephasing times, respectively. Here we have to con- 
sider density matrix elements as the quantities averaged 
over environments degrees of freedom: A = (|C_(i)p|), 
B = lm(C+{t)C*_it)), C = Re{C+it)C*_{t)). The Eqs. 



(|14I15I15|I are similar to well known Bloch-Redfield equa- 
tions in NMRH. 

For spin boson model of coupling of a qubit to thermal 
bath these times have been calculated for weak damping 
Ohmic spectrum by several authors (seei^*i^ and refer- 
ences therein). Here for estimations we take the expres- 
sions for Tr and Td fromi^ 



1 



J{n) coth 



nn 

2fcsT 



Td 2Tr \nJ h ' 



(17) 



(18) 



where dimensionless parameter rj reflects the strength of 
Ohmic dissipation. J{il) is the spectral density of the 
bath fluctuations at the gap frequency. 

The decoherence is caused primarily by coupling of a 
solid state based phase qubit to microscopic degrees of 
freedom in the solid. Fortunately this intrinsic decoher- 
ence has been found to be quite weakii. However, the 
external sources of decoherence are more serious. Here 
we assume that the main source of decoherence is the ex- 
ternal flux noise. For our computer simulations we take 
Tr = 50 ms, Td ^ 200 ns". The evolution of A, B found 
from Eqs. H14I15I16() is shown in Fig. 

As is seen from the Fig. O A decays to 0.5 oscillating 
with Rabi frequency, while B (C) decays to zero. (Note: 
to be rigorous, the stable state solution for A isi^ 



A: 



eq 

Poo 



2(1 + n'„TrTd) 



(19) 



However, as distinct from conventional NMR where 
flj^TrTd ~ 1, for our values of flu, Tr and Td we get 
^"uTrTd « 10'^. For the same reason the stable state os- 
cillations of B and C are quite small). It means the 
qubit density matrix becomes the statistical mixture: 
Poo 1/2, pii 1/2, pio = poi ^ at t ^ oo. 
Therefore, the noise from environment leads to a delo- 
calization: the system which initially is localized in any 
state would be always delocalized at t ^ oo. 
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This property has been first pointed out in connection 
with noninvasive measurements of coherent dynamics in 
quantum-dot systemsiS where one can find an interesting 
discussion of how this behavior is related to the quan- 
tum Zcno effect, and recently has been confirmed for a 
quantum-dot qubit interacting with an environment and 
continuously monitored by a tunnel-junction detector^^. 

As is seen from Eq. [7} , the delocalization leads to the 
vanishing of the current in the qubit loop. 



where 77 = The function F[^{t)] in ^ stands 

for stochastic nature of the measuring process. In ac- 
cordance with von Neumann postulate the outcome of 
a single measurement cannot be predicted deterministi- 
cally. When qubit is in a superposition of two stationary 
states its wave function can be expressed in the flux 
basis as ^{t) ^U'^l + W'^r where 

\U\^ ^ al + Aiai - al) + 2Ca+a-, (27) 



III. THE MEASURING OF FLUX QUBIT WITH 
A TANK CIRCUIT 



Below we consider a continuous measurements of a 
qubit with a classical tank circuit which is weakly cou- 
pled to the qubit via mutual inductance M. First we 
study the ideal case when qubit and tank circuit are 
completely decoupled from their environments, so that 
we may describe qubit -I- tank system by Hamiltonian 
H = Hq + Ht + Hint, where 



(20) 



In (|20|l Ct, and Lt are capacitor and inductor of a 
tank circuit; Q is the charge at the capacitor, $ is the 
magnetic flux trapped by the inductor. The tank-qubit 
interaction is described by Hamiltonian 



Hi, 



(21) 



where qubit current operator Iq is given in ©, A = 
M/Lt. 

The equations for qubit -I- tank system are 

f = 2^cos(.t)S-2A,/cABcI>, (22) 



\W\' 



2C6+6_ 



(28) 



with |C7|2 -f |T^|2 ^ 1. 

The states 5*^, n have equal currents circulating in 
opposite directions so that the outcome of the measure- 
ment (the direction of a current circulation) can be pre- 
dicted only statistically with probability |C/p or 
respectively. Following this reasoning we take -F[C(i)] 
the form: 



mt)] 



\u?-m 
\\u?-m[ 



(29) 



where ^{t) generates random numbers from interval 
[0,1]. Since |f/p < 1 the function F[£,it)] takes two val- 
ues: -|-1, —1. It accounts for unpredictable nature of a 
single measurement. If, for example, in some moment U 
\U{ti)\'^ > 0.5 it is more probable to find the clockwise 
than counterclockwise direction of circulating current at 
this moment of time in a single measurement. The ac- 
tual value of the voltage across the tank at some time 
U is obtained as the average of individual measurements 
over N different realizations of ^ (t) : 



=N 



viu) 



NCi 



(30) 



dB Aei , , 

— = —J-cos(u;t){2A-l)- 

2{n + ^ cos(tJi)) C + 
2A/c^(2eoC + A(2A-l))$, 



(23) 



'^^2{n+'-^cos{ujt))B-2\f^Ic^B>^>, (24) 



(25) 



Therefore, our model accounts for stochastic back ac- 
tion influence of the measuring device (tank circuit) on 
the qubit behavior. In some sense the model resembles 
the probabilistic measurement scheme described in^'^" for 
detection of the electron position in double quantum dot 
by measuring the current through tunnel junction cou- 
pled to quantum double dot qubit. 

Below, the tank circuit parameters are Lt = 50 nH, 
Ct = 200 pF, so that the tank is tuned to 50 MHz. The 
inductive coupling to the qubit M = 12.5 pH that gives 
for the coupling parameter A — 2.5 x lO^''. In addition, 
we take eg = ^ so that 



1 



V2 + I 



A + C 



(31) 



dQ $ 



dt 



Lt 



XvIcFm], 



(26) 



The results of computer simulations of the equation set 
I I26|l are shown in Figs 0] and |S1 
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FIG. 4: Phase loss-free qubit coupled to a loss-free tank cir- 
cuit. Oscillations of A. Deterministic case (a) together with 
one realization (b) are shown. Small scale time oscillations 
correspond to Rabi frequency. 
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FIG. 5: Phase loss-free qubit coupled to a loss-free tank cir- 
cuit. Voltage across the tank. Deterministic case (a) together 
with one realization (b) are shown. 




FIG. 6: Phase loss-free qubit coupled to a dissipative tank 
circuit. The evolution of A exhibits modulation of Rabi oscil- 
lations with lower frequency. Deterministic case (a) together 
with one realization (b) are shown. 




2 4 6 8 10 

FIG. 7: Phase loss-free qubit coupled to a dissipative tank 
circuit. The voltage across the tank exhibits modulation of 
Rabi frequency. Deterministic case (a) together with one re- 
alization (b) are shown. 



At every graph of the figures the results for one realiza- 
tion of random number generator ^ (t) are compared with 
the case when we replaced ™ with determin- 

istic term {2 A — 1 — 2C) / \/2, which means that the tank 
measures the average current {T)) in a qubit loop. As is 
seen from the Fig. 01 A oscillates with Rabi frequency. 
The voltage across tank circuit oscillates also with Rabi 
frequency which is equal to 50 MHz in our case (Fig. 
which is modulated with the lower frequency the value 
of which is about 5 MHz. 

It is worth to note the interesting feature of the result: 
though the system is free from dissipation the voltage 
across the tank is not saturated (the voltage amplitude is 
about 10 nV, Fig.[SJ- Although A oscillates at resonance 
frequency of the tank, the saturation is not reached. A 
simple estimations show that at the tank resonance the 
saturated value of the voltage is about 75 nV. We have 
found the effect is due to the large value of a coupling 
constant A. The simulations show the full saturation is 
reached with A w 10"^, however, then the voltage is quite 
small to be detected. These results are valid exactly for 



loss-free tank circuit, however we may within simulation 
time span (10 /is) extrapolate them for tank with high 
quality factor, say Qt > 1000. 

Now we want to account for the damping in the tank 
circuit. We replace Eq. 1261 with 

^ = -i-i^e-A.*m.)l. (32) 

For the simulations we take tank quality Qt = 100. The 
results of simulations of equation set (|22I231 1241 1251 132(1 
are shown for A in Figs. El and for the voltage in Figs.[7| 

It is worthwhile to note that though the qubit is uncou- 
pled from its own environment, nevertheless, the current 
in a qubit loop and the voltage across the tank decay. 
The quantities B and C which are not shown here os- 
cillate without damping with the frequency of excitation 
with the amplitude equal to 0.5 . 

Finally, we consider the case when the qubit and the 
tank are coupled to their own environments. The corre- 
sponding set of equations are Eqs. (|25I32|) and following 
equations: 
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FIG. 8: Phase dissipative qubit coupled to a loss-free tank 
circuit (1/Qt = 0). The evolution of A. Deterministic case 
(a) together with one realization (b) are shown. 
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FIG. 9: Phase dissipative qubit coupled to a loss-free tank 
circuit {1/Qt — 0). The voltage across tank circuit. Deter- 
ministic case (a) together with one realization (b) are shown. 



(33) 

2 (r! + ^ cos(wt)) C + 
2A/c^(2£oC + A(2A-l))$-^, (34) 

f .2(..£fcos(..))i.-2A,..|^i..-|. 

(35) 

The results of simulations of the equation set 
(|^:^l:^4lj5ll^5l3i>|) are shown in Figs, il El [101 CH and 
Fig. 112 

From Fig. island Fig. ^Iwe see that A is almost unaf- 
fected by the value of Qt- Its decay is defined primarily 
by the shortest time T^. The decay time of the voltage 




FIG. 10: Phase dissipative qubit coupled to a dissipative tank 
circuit [Qt = 100). The evolution of A. One realization 
(lower graph) together with deterministic case (upper graph) 
are shown. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 11: Phase dissipative qubit coupled to a dissipative tank 
circuit (Qt = 100). The voltage across tank circuit. One real- 
ization (lower graph) together with deterministic case (upper 
graph) are shown. 




FIG. 12: Phase dissipative qubit coupled to a dissipative tank 
circuit. The voltage across the tank for deterministic case. 
Quality factors: Qt ~ 1000 (upper graph) and Qt ~ 100 
(lower graph). 
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FIG. 13: Loss-free qubit coupled to a dissipative tank circuit 
(Qr = 100), A = 2.5 X 10"^. The evolution of A. A distance 
between neighbor jumps is equal to Rabi period. Determinis- 
tic case is shown. 




t(^s) 

FIG. 14: Loss-free qubit coupled to a dissipative tank circuit 
(Qt = 100), A = 2.5 X 10"^. The evolution of the voltage 
across the tank for the deterministic case. 



is defined by the value of Qt, since tank circuit decay 
time 2Qt/cot is much longer than the dephasing time 
Td (Figs. 1^ 111! I12|) . It is advantageous from the point 
of experiment, since we can measure Rabi oscillations 
much longer than the dephasing time T^. However, from 
the other hand, the maximum value of the voltage ampH- 
tude is almost independent of Qt (see Fig. I12f) . We have 
shown in^^ that for carefully made electronics the voltage 
noise at the input of preamplifier could be on the order of 
10 nV-20 nV. As is seen from Fig.lllland especially, from 
Fig. El during the first microsecond the signal-to-noise 
ratio is about 0.5. The signal can be recovered with the 
well known in NMR pulsed technique with subsequent 
Fourier processing. However, since here the pulse filling 
frequency is uncoupled from the frequency of the signal 
to be detected, it is not necessary to keep the pulse width 
shorter than the decaying time of the signal. As our re- 
sults show, the pulse duration of about 2 fis is adequate 
for the measurements. 

In conclusion we want to show the effect of qubit evo- 
lution as the coupling between the qubit and the tank is 
increased. We numerically solved the system consisting 
of the loss-free qubit coupled to the dissipative tank cir- 



cuit. The system is described by Eqs. H14I15I16I25() and 

Eq. H32|) . For the simulations we take the coupling pa- 
rameter A = 2.5 X 10^^. The results of simulations are 
shown on Figs. ^1 E| for deterministic case. As is seen 
from the Figs. E] during Rabi period the quantity A be- 
came partially frozen at some level. At the endpoints of 
this period the system tries to escape to another level of 
A. Between the endpoints of Rabi period A oscillates with 
a high frequency which is about 10 GHz in our case. As 
expected, the evolution of B is suppressed approximately 
by a factor of ten below its free evolution amplitude which 
is equal to 0.5. As we show below, the strong coupfing 
completely destroys the phase coherence between qubit 
states, nevertheless the voltage across the tank oscillates 
with Rabi frequency. Its amplitude is considerably in- 
creased and it does not reveal any peculiarities associated 
with the frozen behavior of A (Fis. I14|l . 

IV. QUBIT WAVE FUNCTION 

In conclusion we want to study the effect of a coupling 
between qubit and the tank on the qubit state, in particu- 
lar, on phase coherence between basis states of the qubit. 
It is necessary to note that our measurement is not the 
measurement in the sense of Neumann. We are interested 
only in the voltage amplitude in the tank but not in the 
state of the qubit: we did not solve Schrodinger equations 
for C- = |C_| exp(i(^_), C-^ = |C+| exp(i(^+) but for 
their products A = |C_|^, B = |C_ ||C+| sin((^+ — (y9_), 
C — |C_||C+| cos(i^-|_ — ip^). Nevertheless, we can check 
to what extent the qubit can be described by the wave 
function in case of its interaction with a tank. Evidently, 
free qubit must have definite wave function at any instant 
of time. It means the conservation of phase coherence 
the condition for which can be expressed in terms of our 
quantities as: 

B' + C' _ ^ .... 

^^^^ 

If we switch on the interaction with a tank we may not, 
strictly speaking, consider qubit as having definite wave 
function. However, if the interaction is rather weak the 
qubit wave function could be well defined. We showed 
before that for relatively weak coupling the dissipation 
resulted in quenching A to the 0.5 level (see Figs. El (HI 
I10|l . That means |C+| = |C_| However, as is 

seen from Fig. 1151 the condition of phase coherence is 
still vafid up to A « 10"^. 

As the coupling is increased further the qubit wave 
function is completely destroyed. The quantity A is 
quenched to approximately 0.85 (Fig. I13|l . That is 
|C_| ~ 0.92. It might seem that we have here so called 
Zeno effect- as if qubit state is frozen in its ground state. 
However, in case of a strong coupling it is not correct 
to say about wave function of the qubit alone. This is 
shown in Fig. ll5l where for A > 10^^ the phase coherence 
is seen to be completely lost . 
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FIG. 15: The condition of phase coherence H3()^ vs coupling 
strength A. 



V. CONCLUSION 

In all computer simulations we systematically com- 
pared different realizations with the case when we re- 
placed in with deterministic term {2A — 1 — 
2C)/^/2. We have found that within a decaying time 
all realizations and deterministic case give almost identi- 
cal results. That is why in all corresponding graphs we 
compared deterministic case with the only realization. 
A clear difference appears only at the tails where A is 
close to 0.5 and Cis rather small. This is because at the 
tails the random number generator every moment of time 
changes the sign of the current with a high probability, 
while within a decaying time where A undergoes oscil- 
lations the sign of the current for one half of period of 
oscillations of A is conserved with a high probability. 

Throughout the paper we stress the similarity between 
the qubit-l-tank system and NMR, however, we have to 
be aware of the main difference. In NMR the back ac- 



tion of the tank circuit on the sample under study is 
neglected. It is justified by the fact that the tank is cou- 
pled to macroscopic number of two-level systems (1/2- 
spin particles). The coupling to the individual particle 
is rather small, so that a reasonable signal level at the 
tank is obtained at the expense of enormous number of 
the coupled particles. However, when the tank is cou- 
pled to a single two-level system the account for back 
action is necessary. It leads to the main quantitative dif- 
ference from NMR. In order to keep the noise from the 
tank as small as possible, the quality factor Qt should 
be taken as high as possible from technological point of 
view (in our simulations we take Qt = 100 only in order 
to save the simulation time). However, the signal ampli- 
tude weakly depends on Qt being at best at the level of 
noise. Nevertheless, it is not difficult to recover the signal 
with the aid of the methods of signal processing which 
are used in high resolution NMR. Therefore, the results 
of our simulations clearly show that we can detect Raby 
oscillations of the voltage across tank circuit coupled to 
the qubit with the pulsed Fourier technique which is well 
known in NMR. 
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